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1.0  INTRODUCTION 

Natural  water  droplets  and  ice  crystals  in  the  atmosphere  present  a  formidable  hazard  to 
ballistic  vehicles  during  various  phases  of  their  trajectory.  A  number  of  studies  have  been 
made  to  determine  the  effect  of  ice  crystals  and  water  droplets  on  the  performance  of 
ballistic  reentry  vehicles  (RV)  (Refs.  1  and  2).  The  gun  range  has  been  used  to  determine  the 
effect  of  snow  and  ice  crystal  impact  on  reentry  vehicle  nosetips  for  extremely  short  flight 
times  (Ref.  3).  At  present,  there  is  no  facility  that  can  simulate  both  the  duration  and 
intensity  of  the  erosive  environment  induced  by  natural  water  droplets  during  the  reentry 
phase  of  ballistic  flight  trajectories.  The  Arnold  Engineering  Development  Center  (AEDC) 
Dust  Erosion  Tunnel  (DET)  has  been  used  successfully  to  simulate  the  environment  created 
by  atmospheric  dust  particles  during  the  ascent  phase  of  a  ballistic  vehicle’s  flight.  In  the 
AEDC  High  Enthalpy  Arc  Tunnel  (HEAT)  test  facility  (H-l)*,  the  kinetic  energy  flux  on 
transpiration-cooled  nosetips  attributable  to  water  droplet  impacts  has  been  simulated  with 
graphite  particles  having  a  specific  gravity  of  1.7.  Extending  the  simulation  capability  to 
include  acceleration  of  water  droplets  would  provide  a  more  accurate  simulation  of  the 
environment  encountered  by  an  RV  nosetip  in  adverse  weather  during  its  ascent  phase  in  the 
DET  and  its  terminal  phase  in  the  HEAT  facility.  The  survivability  of  droplets  in  these  high- 
temperature  facilities  is,  of  course,  the  critical  feasibility  question. 

In  Fiscal  Year  1977,  an  exploratory  test  was  conducted  by  H.  F.  Lewis  at  AEDC  in  the 
DET  in  which  liquid  water  (up  to  15  percent  of  total  mass  flow)  was  injected  through  the 
existing  particle  injection  system.  The  facility  was  operated  at  low  enthalpy  (600  Btu/lbm  to 
encourage  droplet  survival.  No  particle  diagnostics  were  installed  to  provide  direct  evidence 
of  droplet  survival;  however,  strong  indirect  evidence  was  obtained  from  stagnation  point 
heat-transfer  measurements  which  could  be  interpreted  in  terms  of  the  latent  heat  of 
vaporization  of  droplets  forming  a  heat  sink  at  the  model  surface.  In  addition,  some 
evidence  of  model  surface  erosion  and/or  accelerated  chemical  activity  was  observed.  The 
current  study  w'as  undertaken  to  determine  whether  droplet  survival  in  the  DET  tests  could 
be  analytically  verified  and,  more  importantly,  to  determine  whether  appropriate 
modifications  would  produce  a  useful  droplet  erosive  field  in  the  DET  and  the  HEAT 
facility. 


2.0  DROPLET  BEHAVIOR 

This  report  is  concerned  with  understanding  and  predicting  the  behavior  of  accelerating 
water  droplets  in  supersonic  nozzles  at  operating  conditions  usually  associated  with 
complete  test  model  destruction.  For  a  water  droplet  to  survive  such  conditions,  it  must  meet 
and  overcome  three  specific  hazards: 


•The  AEDC  High  Enthalpy  Arc  Tunnel  (HEAT)  test  facility  (H-l)  is  hereinafter  referred  to  as  the  HEAT  facility. 
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1.  The  droplet  or  its  precursor  must  be  injected  into  the  nozzle  in  a  manner 
consistent  with  the  formation  of  the  droplet  and  the  basic  flow  process  in  the 
nozzle. 

2.  The  droplet  must  resist  breakup  by  the  aerodynamic  forces  that  produce  its 
acceleration. 

3.  Enough  of  the  initial  droplet  must  survive  evaporation  to  produce  a  meaningful 
test  environment. 

The  design  goal  for  a  droplet  acceleration  nozzle  is  to  accelerate  water  droplets  of  from 
500-  to  1 ,000-^  diameters  to  velocities  that  are  associated  with  ascent  and  reentry 
trajectories.  The  operating  parameters  for  the  DET  and  HEAT  are  presented  in  Table  1. 

Since  there  is  no  continuous  water  droplet  acceleration  facility  presently  available, 
attainment  of  even  a  fraction  of  the  velocities  indicated  in  Table  1  will  be  of  some  interest.  It 
is  expected,  however,  that  droplet  velocities  will  have  to  exceed  9,000  fps  to  be  of  lasting 
interest.  Also,  although  droplet  sizes  of  500  fi  are  desired,  more  interest  would  result  if 
smaller  drops  could  be  accelerated  to  high  velocities. 


Table  1.  DET  and  HEAT  Facility  Operating  Parameters. 


DET 

HEAT 

Throat  Diameter  (in.) 

0.5625 

0.375  to  0.700 

Exit  Diameter  (in.) 

4.36,  8.32,  15.34 

0.850  to  1 .600 

Mach  Number 

5.9,  7.7,  9.5 

1.8  to  3.5 

PQ  (atm) 

20  to  68 

20  to  120 

hQB  (Btu/lmb) 

400  to  1 , 800 

400  to  4,000 

Pe  (atm) 

0.044  to  0.0006 

0 . 50  Min . 

u  (fps) 

gas  '  r  ' 

4,300  to  9,200 

2,800  to  11,100 
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2.1  DROPLET  FORMATION 

Since  droplet  formation  has  been  of  significant  industrial  importance  for  many  years, 
there  has  been  much  research  done  concerning  the  breakup  of  liquid  streams  into  drops. 
Although  this  study  is  restricted  to  the  study  of  droplet  injectors  for  the  DET  and  HEAT 
nozzles,  a  short  discussion  of  the  full  range  of  liquid  jet  breakup  behavior  will  illuminate 
some  of  the  problems  inherent  in  droplet  injector  design. 

In  Ref.  4  and  others,  the  jet  breakup  behavior  is  illustrated  by  showing  a  simplified 
model  of  an  ideal,  smooth,  cylindrical  injection  tube  discharging  liquid  into  stagnant  air  at 
various  liquid  How  rates.  The  resulting  flow  patterns  for  an  increasing  liquid  flow  rate  are 
illustrated  in  Figs,  la  through  e. 
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Figure  1.  Observed  liquid  jet  breakup  regimes. 
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The  flow  pattern  sketched  in  Fig.  la  represents  a  low  flow  regime  in  which  the  capillary 
forces  completely  dominate  the  flow.  A  drop  remains  attached  to  the  injector  until  it  is 
dislodged  by  the  weight  of  the  fluid,  stray  air  currents,  vibration,  or  some  other  force.  The 
size  of  the  drop  is  controlled  by  the  liquid  surface  tension  and  injector  geometry.  Many 
experimenters  have  reported  that  the  drop  size  is  roughly  twice  the  injector  inside  diameter. 

In  Fig.  lb,  the  liquid  flow  rate  has  been  increased  until  a  solid,  clear,  undisturbed  jet  is 
formed.  The  clarity  of  the  jet  is  significant  in  that  there  are  no  surface  waves  on  the  jet 
surface;  this  implies  that  the  aerodynamic  forces  are  minimal.  The  jet  surface  oscillates 
because  of  the  interaction  between  surface  tension  and  inertial  forces;  these  surface 
oscillations  build  to  the  point  at  which  the  jet  is  pinched  off,  forming  long  cylinders 
downstream  of  the  injection  point.  These  cylinders,  in  turn,  form  droplets  that  are 
approximately  twice  the  jet  diameter.  This  regime  is  called  the  Rayleigh  regime  and  has  been 
studied  intensively. 

At  higher  jet  velocities,  the  aerodynamic  forces  induced  by  the  motion  of  the  jet  become 
significant.  The  aerodynamic  forces  amplify  and  interact  with  the  surface  oscillations  to 
increase  their  amplitude  and  produce  jet  breakup  nearer  the  injector  face,  as  is  shown  in  Fig. 
Ic.  The  resulting  droplets  have  the  same  diameter  as  the  jet.  This  regime  was  called  the  “first 
wind-induced  regime”  by  Weber,  whose  work  is  discussed  in  Refs.  5  and  6. 

The  second  wind-induced  regime,  which  occurs  at  even  higher  jet  velocities,  is 
characterized  by  a  relatively  short  cylindrical  jet  that  abruptly  breaks  up  into  small  drops, 
producing  a  fan-  or  cone-like  expansion,  as  shown  in  Fig.  Id.  In  this  regime  the  surface 
aerodynamic  forces  are  so  much  greater  than  the  surface  tension  forces  that  the 
characteristic  drop  size,  which  is  considerably  smaller  than  the  jet  diameter,  is  only  slightly 
dependent  upon  jet  diameter. 

The  final  regime,  shown  in  Fig.  le,  is  called  the  atomization  regime  because  the  jet  is 
completely  broken  up  into  very  small  drops  (i.e. ,  atomized)  at  the  injector  exit.  The  jet  has 
no  cylindrical  section  and  the  cone-like  expansion  begins  at  the  injector  exit.  Some 
experimental  evidence  suggests  that  the  breakup  process  actually  begins  inside  the  injector 
and  that  it  is  associated  with  turbulence  production  in  the  injector. 

Much  empirical  and  analytical  work  has  been  done  in  attempts  to  describe  the  jet 
breakup  in  the  atomization  regime;  important  applications  to  such  devices  as  paint  sprayers 
and  rocket  propellant  injectors  have  helped  to  prompt  this  work.  There  is  a  lack  of 
analytical  definition  of  the  boundaries  of  the  various  breakup  regimes,  probably  because  the 
many  secondary  variables  that  affect  the  breakup  process  are  uncontrolled,  thus 
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complicating  the  precise  definition  of  boundaries.  The  Rayleigh  and  first  wind-induced 
breakup  regimes  are  of  primary  importance  to  this  study;  it  is  unfortunate  that  the  other, 
more  stable,  breakup  regimes  produce  large  quantities  of  very  small  droplets  which  tend  to 
evaporate  quickly  in  a  nozzle. 

The  two  critical  questions  of  droplet  behavior  that  must  be  answered  are 

1.  Where  does  the  liquid  jet  breakup  occur? 

2.  What  is  the  diameter  of  the  resulting  drops? 

Calculation  of  a  meaningful  jet  breakup  time  in  either  the  DET  or  the  HEAT  facility  is 
complicated  by  the  fact  that  the  injected  jet  is  subjected  to  an  extreme  acoustic  environment 
and  a  highly  turbulent  flow.  Measurements  of  the  environment  inside  the  injector  chamber 
have  not  been  made,  but  pitot  pressure  measurements  in  the  arc  heater  exit  flow  typically 
have  fluctuations  that  are  2.5  percent  of  the  mean  value.  The  quantitative  effect  of  the 
droplet  injection  chamber  environment  is  also  unknown.  In  view  of  the  uncertain  nature  of 
the  environment  and  its  effect  on  the  jet,  the  most  conservative  assumption  to  make  is  that 
the  jet  breaks  up  into  droplets  at  the  injector  exit,  although  it  could  be  possible  that  the  jet 
could  break  up  further  downstream,  producing  larger  droplets  at  the  exit. 

Many  different  droplet-size  prediction  models  exist.  The  experimental  work  of  Ingebo 
(Ref.  7)  is  similar  to  the  earlier  DET  test  in  terms  of  type  of  injection,  injection  liquid, 
flow-field  variables,  and  gas  and  liquid  velocities,  and  is  therefore  most  applicable  to  this 
study.  Ingebo  has  correlated  his  results  for  accelerating  water  and  ethanol  drops  in  still  and 
accelerating  air,  nitrogen,  and  helium  streams  as 

B0-33  h  Rep01  Re0  5  B'°  07 

o  l  g  O 

(l.25  X  lO'6  +  0.33  X  106  A°  5  13;°  67)  J  0) 


O  R 


D 


=  0.64 


'  n  \0  33 

UOR  \ 


Max 


[o.044  *WeEWeg02B-°-07 
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P£  dor»5 


n 
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Ke  -  Weber  Number  — 
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P,  Vr  DOH 


(3) 
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Pg  D0  R  a 

A  =  Acceleration  to  surface  tension  ratio  =  - 

C 


(4) 


The  remaining  terms  and  subscripts  are  explained  in  the  Nomenclature. 
The  drop  or  jet  velocity  issuing  from  an  injector  can  be  predicted  by 


where  <5P  is  the  injector  pressure  drop  and  Cq  is  the  empirical  discharge  coefficient,  usually 
between  0.85  and  1 .0. 

A  plot  of  Eq.  (5)  is  shown  in  Fig.  2  for  Cq  =  1.0.  It  should  be  noted  that  in  all  of  the 
literature  searched  concerning  “high-speed  jets,”  injector  5P  is  limited  to  about  3,000  psi. 
No  data  and,  hence,  no  correlations  were  found  for  what  may  be  termed  ultrahigh  jet 
velocities  on  the  order  of  2,000  to  3,000  fps. 


0  20,000  40,000  60,000  80.000  100,000 

Injector  Pressure  Drop,  ap,  psi 


Figure  2.  Theoretical  water  droplet  velocity  as  a 
function  of  injector  pressure  drop. 
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Discussion  of  the  practical  problems  associated  with  achieving  high  jet  velocities  will  be 
deferred  until  Section  4.3,  but  it  is  appropriate  to  discuss  the  jet  breakup  problem  here. 
Some  investigators  believe  that  atomization  of  the  jet  occurs  because  of  the  level  of 
turbulence  in  the  injector  and  not  because  of  the  jet  velocity  relative  to  the  airstream.  Figure 
3  presents  the  maximum  droplet  diameter,  DMax,  predicted  by  Eq.  (1)  as  a  function  of  the 
injection  orifice  diameter.  Dor,  for  two  droplet  injection  chamber  conditions. 


Jet  Injection  Velocity,  fps 

Figure  3.  Predicted  maximum  droplet  diameter  as  a 
function  of  jet  injection  velocity. 

2.2  DROPLET  BREAKUP 


The  most  restrictive  feature  of  the  droplet  behavior  is  the  tendency  for  the  droplet  to 
deform  and  disintegrate  when  the  aerodynamic  loading  exceeds  a  critical  loading.  The  ratio 
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of  droplet  distorting  aerodynamic  forces  to  the  droplet  cohesive  forces  is  defined  as  the 
Weber  number: 


We 


P ,  <U~V)2  Dd 


(6) 


A  number  of  experimental  stidoes  of  these  phenomena  have  been  made  (Refs.  9,  10,  and  1 1). 
The  behavior  of  the  drop  as  described  in  Ref.  4  is  categorized  by  the  five  modes  shown  in 
Fig.  4. 
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Ref.  9 


Figure  4.  Observed  droplet  breakup  regimes. 


The  critical  Weber  number  is  defined  as  the  upper  boundary  of  the  first  region,  that  of 
steady-state  equilibrium.  Values  vary  from  4  to  10  for  the  critical  Weber  number  used  in  the 
studies.  The  value  of  8,  taken  from  Ref.  9,  is  used  in  the  analytical  portion  of  this  study. 
Between  the  critical  Weber  number  and  a  Weber  number  of  20,  the  droplet  is  observed  to 
vibrate  and  eventually  to  break  up  into  two  or  three  equal-sized  droplets.  Between  Weber 
numbers  20  and  70,  the  droplet  deforms  into  a  bag  with  a  heavy  rim  before  breakup. 
Between  Weber  numbers  70  and  100,  the  drop  deforms  as  a  bag  plus  a  centerbody  before 
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breakup.  Between  numbers  100  and  20,000,  the  drop  flattens  and  the  surface  is  blown  (i.e., 
stripped)  away.  Above  a  Weber  number  of  20,000,  the  drop  is  said  to  enter  the  catastrophic 
mode  of  breakup  in  which  it  is  torn  apart  virtually  instantaneously.  Two  points  should  be 
apparent  from  Fig.  4:  (1)  Drops  deform  before  breakup,  and  {2)  droplet  breakup  is  a 
transient  process  below  a  Weber  number  of  20,000. 


Reinecke  and  Waldman  performed  a  number  of  droplet  experiments  with  a  shock  tube; 
they  observed  the  behavior  of  droplets  with  flash  x-ray  and  shadowgraph  devices.  They 
correlated  (Ref.  9)  the  time  required  for  a  droplet  to  break  up  after  passing  through  a  normal 
shock  wave  as 

220*  -  e)'* 

T,  =  - — - For  8<W  <1,700 

o  i,  e  _  ’ 

QMai 1 


where 


0.013(W  —  8)-'* 

Tb  =  - - -  For  1,700  <  <  2,700 

Q|U  ax 


35.0(W  -  8r-,« 

T.  =  — - - -r, -  For  7,700  <  W 

Qm 


Qm£ 


h  x  on  D  ro p 

%c 


=  0.78  + 


1.47 


1  0  *■  2.1  M34 


(7) 


(8) 


T 


b 


Lb  Vr  Pe 

Dc.  Pi 


(9) 


In  these  experiments  the  flow  behind  the  shock  is  steady,  and  the  droplet  is  accelerating 
toward  equilibrium  with  the  flow.  In  the  droplet  acceleration  nozzles  under  study  here,  both 
the  droplet  and  the  gas  are  accelerating.  Equations  (6),  (7),  and  (8)  are  used  as  the  transient 
droplet  breakup  criteria  by  requiring  that 

1Th 

—  =  1 

Tb  (10) 

before  the  droplet  can  be  said  to  have  disintegrated. 
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The  droplet  dilation  ratio  observed  in  Ref.  9  can  be  summarized  as  follows: 


Dj  2 

-  =  1  +-  1.51  1.75  Max.  observed 


For  0  <  W  <  2,000  (11) 


—  =  1  +  2.2T  4.8  Max,  observed  0  <  T  <  1.8  For  2,000  <  ffe  <  3,200  (12) 


-  1.2031  3.5  Max.  observed 


For  3,200  <  W 


(13) 


The  dilation  behavior  in  the  intermediate  Weber  number  region  was  odd,  and  the  scatter 
in  all  regions  was  large.  The  relations 


D. 

—  =  1+  1.5T2 
D, 

d  i 


(14) 


3) 

Max 


2.5 


(15) 


are  used  in  the  study  for  all  Weber  numbers. 


2.3  DROPLET  DRAG  COEFFICIENT 


According  to  Ref.  12,  the  net  force  on  a  droplet  is  comprised  of  several  components: 


F 


F  _ 


Ba 


F 

am 


+ 


Hod  v 


(16) 


where 


Ft  =  Viscous  shear  forces 

Fp  =  Surface  pressure  forces 

Fam  =  Force  attributable  to  mass  entrained 
in  droplet  boundary  layer 

F0a  =  Bassett  force  attributable  to  unsteady 
motion  of  flow  around  the  droplet 
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Fbody  =  Body  force  attributable  to  gravity  and 
to  electrostatic  and  magnetic  fields. 

The  forces  attributable  to  gravity,  electromagnetic  fields,  apparent  mass,  and  unsteady 
fluid  motion  can  be  shown  to  be  small  for  droplet  and  nozzle  flow  conditions  under 
consideration.  |See  Soo  (Ref.  12)  and  Rudinger  (Ref.  13)  for  details.] 

The  electrostatic  forces  can  also  be  neglected  in  this  study,  but  under  some  conditions, 
the  electrostatic  force  either  advertently  or  inadvertently  applied  can  be  significant  (Ref.  14). 
The  remaining  forces  are  grouped  together  and  are  described  by  the  droplet  drag  coefficient 
to  which  appropriate  corrections  are  made  to  account  for  the  effects  of  nozzle  pressure 
gradient  and  droplet  volume  fraction. 

Comparison  of  droplet  drag  coefficient  measurements  is  complicated  by  the  tendency  of 
the  droplets’  surfaces  to  distort,  slip,  and  vaporize.  A  comparison  of  various  data  to  the 
standard  incompressible  sphere  drag  curve  is  shown  in  Fig.  5,  taken  from  Ref.  15.  The  data 
of  Hughes  and  Gilliland  (Ref.  16)  are  for  freely  falling  drops  of  various  diameters.  The  data 
of  Ingebo  (Ref.  17)  are  for  accelerating  droplets  of  various  liquids.  The  disagreement  or 
scatter  of  these  data  (i.e.,  from  0.2  to  1.0),  at  Reynolds  numbers  between  400  and  500 
indicates  that  droplet  calculations  in  this  range  are  approximate  at  best. 


Figure  5.  Drag  coefficient  of  spheres  and  drops. 
Ref.  16. 


The  curve  fit  to  a  standard  drag  curve  presented  in  Ref.  15 


Cd  =  —  For  Red<0.01 
Re, 


(17) 
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r  ,,  2 

C,  =  exp  3.271  -0.8893  (£n  Hcd)  +  0.03  117  ((in  Red) 

L  (18) 

+  0.00 1 143  (£n  Re(j)3l  For  0.1  <  Hed  <  20,000 
Cd  =  0.4569  For  Re<[  >  20,000  (19) 


is  used  in  this  study  with  appropriate  correction  factors.  When  one  derives  optimum  nozzle 
contours,  the  use  of  a  single  relationship  for  drag  coefficient  between  Red  =  0.01  and 
Red  =  20,000  precludes  any  problems  attributable  to  drag  curve  mismatch.  Mach  number 
correction  to  the  standard  drag  curve  is  not  necessary  since  the  relative  Mach  number  of 
viable  drops  will  be  low'  in  areas  in  which  small  changes  in  the  drag  coefficient  are 
significant. 


The  nominal  value  of  droplet  Cj  must  be  increased  if  the  droplet  spacing  approaches  a 
value  of  several  droplet  diameters.  Tam  (Ref.  18)  estimates  that  the  drag  coefficient  of  a 
sphere  in  incompressible  flow  will  be  multiplied  by  a  factor  as  follows: 


cd  _  It  3f.  3(.8f  -  3e2)  * 
C,  “  (2  —  3f)z 

Sink  e  s 


(20) 


where  e  is  the  fraction  of  space  occupied  by  the  particle  volume.  The  liquid  volume  fraction 
is  less  than  0.2  percent  for  the  conditions  considered  in  the  study;  consequently,  it  is  not 
necessary  to  make  the  drag  coefficient  correction  indicated  by  Eq.  (20). 


However,  an  additional  correction  is  required  to  account  for  the  effect  of  the  nozzle 
pressure  gradient.  Since  the  force  on  a  droplet  attributable  to  pressure  gradient  is 
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(21) 
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This  correction  is  usually  neglected;  however,  the  small  throat  diameter  and  radius  of 
curvature  and  the  high  reservoir  pressure  of  the  DET  and  HEAT  facility  nozzles  create  an 
unusually  large  pressure  gradient  that  can  increase  the  droplet  drag  coefficient  by  100 
percent.  Because  of  the  nature  of  this  force  (i.e.,  its  similarity  to  the  buoyant  force  which 
accelerates  vapor  bubbles),  it  probably  does  not  increase  the  distortion  of  the  droplet  and, 
hence,  is  highly  advantageous.  Paradoxically,  contouring  the  nozzle  to  produce  a  gradual 
droplet  acceleration,  inhibiting  droplet  breakup  by  aerodynamic  forces,  reduces  this 
favorable  force  to  an  insignificant  level. 

2.4  DROPLET  EVAPORATION 

In  general,  accurate  computation  of  droplet  evaporation  requires  full  consideration  of 
heat  and  mass  transfer  through  the  binary  boundary  layer  of  the  evaporating  drop.  The 
parameters  that  best  describe  this  process  are  the  Nusselt  number  (Nu)  and  the  diffusion 
Nussell  number  (Nu').  These  parameters  have  been  correlated  by  many  investigators  who 
agree  with  Ref.  19  that 


Nu  =  2.0  +  0.60  Re^  Pr*7’ 


and 


Nu'  =  2.0  +  0.60  Re‘J  Sc1'3 


where 


and 


Nu 


D<tht 
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R  T 

Nu'  =  D  ,  — n — — 
d  M 
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g 


(23) 


(24) 


(25) 


(26) 


It  should  be  noted  here  that  Eqs.  (23)  through  (26)  reflect  the  value  of  the  Nusselt  numbers 
observed  in  many  different  tests  and,  as  such,  cannot  predict  the  effects  of  some  aspects  of 
the  flow  that  are  unique  to  the  arc  jet  facilities,  such  as  free-stream  turbulence.  The  heat- 
transfer  and  mass-transfer  rates  to  the  drop  may  be  calculated  from 
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and 


[ 


A  +  Cpv{Tg-Tds)] 


Kg  p  & 


(27) 


(28) 


Equations  (27)  and  (28)  are  also  taken  from  Ref.  19,  although  Eq.  (28)  has  been  modified 
by  the  inclusion  of  the  vapor  pressure  of  the  surrounding  gas  flow,  pVoo.  The  mass  flow  in 
Eq.  (27)  is  the  rate  at  which  the  droplet  loses  mass,  whereas  the  mass  flow  in  Eq.  (28)  is  the 
rate  at  which  the  droplet  mass  is  passed  through  the  droplet’s  boundary  layer  by  diffusion. 
These  rates  may  be  equated  under  most  conditions.  The  energy  balance  of  Eq.  (27)  includes 
a  droplet  radial  heat  conduction  term,  Kf  (3T/dr).  This  usually  insignificant  term  plays  an 
important  role  in  the  change  of  phase  of  a  drop  when  there  is  no  heat  transfer  or  diffusion 
through  the  drop’s  boundary  layer,  such  as  when  a  drop  is  in  both  velocity  and  thermal 
equilibrium  with  a  gas  whose  pressure  is  decreased  to  a  level  below  the  equilbrium  saturation 
pressure.  In  such  a  case,  the  droplet  conditions  necessary  for  evaporation  (i.e., 
nonequilibrium),  are  created  as  the  pressure  wave  associated  with  the  pressure  change  passes 
through  the  drop.  Evaporation  or  boiling  first  occurs  at  a  nucleation  site  or  phase  boundary 
as  the  driving  potential  is  increased.  In  this  case,  the  evaporation  occurs  immediately  at  the 
drop’s  surface  without  any  accompanying  mass  diffusion  through  the  boundary  layer.  The 
evaporation  is  limited  only  by  the  rate  at  which  the  latent  heat  of  vaporization  can  be 
supplied  by  conduction  from  the  interior  of  the  drop.  That  is,  Eq.  (26)  degenerates  to 


m 

Sj  A  dr 


(28) 


Very  large  rates  of  droplet  temperature  change  are  generated  by  seemingly  insignificant 
droplet  temperature  gradients  for  small  drops.  For  example,  the  temperature-time  derivative 
for  a  25-fi  drop  with  a  1  °R  linear  temperature  differential  between  the  surface  and  the  center 
of  the  drop  is 

dj 

dt 


6  kE  dT 

cnw  Dd  dr 


—  =  2,756  °R/ sec 
di 


(29) 


In  such  cases  it  is  plausible  to  assume  that  a  state  of  equilibrium  exists.  This  assumption 
removes  the  need  for  the  radial  conduction  term  in  the  droplet  energy  equation  [Eq.  (27)] 
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which  is  required  only  for  this  case.  The  assumption  can  be  used  to  delete  the  diffusion 
equation  [Eq.  (27)]  and  replace  it  with  the  specification  of  droplet  surface  pressure: 

pvs  =  p,  which  implies  Tds  =  Tvg  =  f(p)  (30) 

where  f  is  the  relationship  of  vapor  temperature  to  vapor  pressure,  with  which  the  droplet 
energy  equation  can  be  solved  for  the  evaporation  rate  as  a  function  of  pressure  and  heat- 
transfer  rate. 


3.0  ANALYSIS  OF  DROPLET  ACCELERATION 

A  review  of  the  earlier  DET  test  data  indicated  that  a  droplet  erosive  field  might  be 
obtained  in  the  DET;  however,  subsequent  analysis  using  the  one-dimensional  particle 
acceleration  nozzle  program  of  Ref.  8,  which  was  modified  to  calculate  droplet  evaporation 
and  breakup,  did  not  support  this  contention.  Further  review  of  the  test  data  and  the 
describing  equations  indicated  that  including  the  momentum  and  energy  coupling  term 
between  droplet  phase  and  nozzle  gas  phase  would  improve  the  analysis. 

3.1  GOVERNING  EQUATIONS 

The  one-dimensional  channel  flow  of  a  noncondensable  gas  with  concurrent  flow  of  a 
liquid  and  its  vapor  can  be  described  by  three  global  conservation  equations,  two  coupling 
relations  for  the  gaseous  and  liquid  phases  momentum  and  energy,  an  equation  describing 
the  rate  of  phase  change,  and  three  equations  of  state.  Auxiliary  relations  are  required  to 
determine  the  global  equation  of  state  in  terms  of  the  component  parts  and  the  gas  phase 
transport  properties  in  terms  of  the  noncondensable  gas  properties,  the  vapor  properties, 
and  the  vapor  fraction,  rha/nig.  The  equations  to  be  integrated  are  as  follows: 

Conservation  of  Global  Mass  [Eq.  (A-2)] 

1  du  ]dh  2^N  1  dp  IdR 
u  dx  h  dx  dx  p  dx  ^  dx 

Conservation  of  Global  Momentum  [Eq.  ( A-4)] 

du  +  *K_  dp+  (u  ~  V 

dx  m  dx  m  dx  m  D, 

8  E  «  d 


dD  , 


D  ,  m  dx 

d  g 


(31) 


dx 


2u  C. 


(32) 
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Conservation  of  Energy  [Eq.  (A-5)J 


du  dh  «"d  dVd  „  dTd\  3  L  dDd  /„  m  V*  „2\  AC. 

Vd17  +  CpJ7r)  +  ~5'17(ci-4T'i+  T~h“7/  *  --r(h»"h* 

£  2  d  4  ^  uv 


(33) 


Particle  Velocity  (from  Newton’s  Second  Law,  applied  to  a  drop)  [Eq.  (15)  of  Ref.  19] 

dV 


d  =  3  C.ptu-V,)  u-V.j 
dx  4 


Dd  Vd  Pd 


(34) 


Particle  Temperature  (from  a  balance  of  aerodynamic  heat  transfer,  internal  energy  of  a 
droplet,  and  the  latent  heat  and  sensible  heat  of  an  evaporating  fraction  of  the  drop)  [Eq. 
(37)  of  Ref.  15] 


where 


(35) 


and 


Hate  of  Heat  Transfer  to  a  Drop 

+  C  (j  -  T  > 

p,  v  °d  v*p) 


(36) 

(37) 


dl,ld 

dx 


—  n 


The  term  (1/v)  (dR/dx)  is  evaluated 
vanishes  for  ideal -gas  calculations. 
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tn 


P  ~  P- 


P  ~  P„ 


(38) 


numerically  from  the  previous  calculation.  This 


The  integration  of  the  equations  is  accomplished  by  estimating  the  downstream  pressure 
and,  hence,  the  pressure  gradient  required  to  obtain  some  desired  downstream  condition, 
such  as  in  matching  the  flow  diameter  to  fit  the  nozzle  contour: 
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dp  =  pu2  J_  (40) 

dx  t_M2  Ak  dx 

Equation  (39)  is  derived  from  isentropic  flow  theory.  This  equation  is  used  successfully 
to  predict  the  value  of  dp/dx  for  nonisentropic  flow  because,  in  practice,  the  5A  term 
includes  the  correction  for  any  previous  errors  and  simplifying  assumptions  as  well  as  for  the 
change  in  true  nozzle  area  in  the  interval  5X. 

This  equation  represents  a  change  from  the  procedure  used  in  Ref.  8  and  has  two  distinct 
advantages: 

1.  In  Eq.  (40)  the  problem  of  dealing  with  the  implicit  sonic  flow  singularity  in  the 
system  of  ordinary  differential  equations  is  replaced  with  the  problem  of  dealing 
with  an  explicit  singularity  whose  behavior  can  be  monitored.  As  a  practical 
matter,  limiting  the  pressure  gradient  to  small  negative  values  in  the  region  near 
Mach  1  eliminates  the  mathematical  problem.  The  gas-dynamic  problem  of 
matching  the  sonic  point  to  the  throat  is  handled,  as  always,  by  iterating  on  inlet 
mass  flow  for  constant  inlet  total  pressure  and  total  enthalpy. 

2.  Equation  (37)  can  be  replaced  with  other  algebraic  equations  to  optimize  the 
nozzle  contour  for  any  conceivable  criteria.  (See  Section  4.1  for  specific  details). 

The  skin-friction  coefficient  for  a  smooth-walled  pipe  is  given  by  the  relation 

C{  =  0,0014  -  0.125  x  Re"0,32  (41) 

taken  from  Ref.  26.  The  nozzle  wall  heat-transfer  rate  is  calculated  with  use  of  Reynolds’ 
analogy  between  heat  and  momentum  transport  as 

Ch  =  0.53  Cf  (42> 

The  equation  of  state  for  air 

p  =  p(p,  h)  (43) 
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is  determined  from  curve  fits  to  the  data  of  Brahinsky  and  Neel  (Ref.  20).  Similarly,  the 
equation  of  state  of  water  vapor 


*\ap  =  (44) 

and  the  definition  of  the  liquid-vapor  interface 

Tvp  -  T<P>  (45) 

pva))  -  Prn  m 

=  h(p)  (47) 

v  ap 

are  determined  by  curve  fits  to  the  data  of  Ref.  21 .  The  surface  tension  of  liquid  wate  in  air 

°  =  ff(T)  (48) 

is  determined  from  curve  fits  to  the  data  of  Refs.  23,  24,  and  25.  The  diffusion  coefficient 
for  water  vapor  in  air 


al)  ®  a  b  rnax’  P^ 


was  determined  from  the  data  of  Ref.  23. 


(49) 


Equations  (31)  through  (49)  are  written  in  terms  of  the  average  gas  and  particle 
properties  and,  as  a  result,  are  as  accurate  as  the  values  of  Cf,  Ch,  and  Cd  used.  With  this  one 
limitation,  there  is  no  reason  why  this  solution  should  not  be  applicable  to  calculations 
beyond  the  merge  point  of  the  turbulent  boundary  layer  of  long  slender  nozzles.  The 
previous  analysis,  presented  in  Ref.  8,  used  a  friction  coefficient  (Cf)  that  was  calculated 
with  the  implicit  assumption  that  the  boundary  layer  was  not  merged;  hence,  this  analysis 
was  optimistic.  The  present  analysis  does  not  make  this  assumption;  however,  in  recognition 
of  the  fact  that  the  centerline  values  of  the  flow  variables  are  different  from  the  average 
values,  the  centerline  flow  conditions  are  estimated  with  the  assumption  that  isentropic  flow 
exists  between  the  nozzle  inlet  and  the  nozzle  exit  along  the  nozzle  centerline.  These 
calculations  are  made  to  determine  to  what  extent  the  nozzle  centerline  conditions  differ 
from  the  average  conditions.  The  resulting  centerline  calculations  are  accurate  up  to  the 
merge  point.  The  centerline  calculations  beyond  the  boundary-layer  merge  point  represent  a 
theoretical  limit  that  cannot  be  reached. 
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3.2  DUST  EROSION  TUNNEL  PREDICTIONS 

The  equations  described  in  Section  3.1  were  programmed  for  an  IBM  370/150  computer 
and  were  used  to  predict  the  droplet  trajectories  in  the  DET  nozzle  for  the  conditions  of  the 
earlier  test: 


p„  =  506  psia 
h„  -  487  Btu/tbm 


rriHjo  =  0.347  Ibm/sec 


The  results  of  the  analysis  suggest  that  small  droplets  could  have  survived  in  the  DET 
nozzle  if  the  droplet  model  described  in  Section  2.0  is  accurate.  The  calculations  of  droplet 
evaporation  predict  that  all  droplets  having  an  initial  diameter  (D.)  greater  than  10  n  will 
survive  evaporation,  as  shown  in  Fig.  6.  Conversely,  the  calculations  of  droplet  breakup 
based  on  the  approximate  model  of  Section  2.2  predict  that  droplets  with  an  initial  diameter 
of  15  g  will  exceed  a  Weber  number  of  8  and,  hence,  begin  to  break  up  0.60  in.  upstream  of 
the  DET  throat,  as  shown  in  Fig.  7;  the  droplets  will  have  completed  the  breakup  process  at 


Figure  6.  Calculated  droplet  diameter  decrease  for  DET. 
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0.15  in.  upstream  ol  the  throat,  as  shown  in  Fig.  8.  The  entire  breakup  process  occurs  in  0.45 
in.  of  droplet  travel,  or  approximately  762  initial,  droplet  diameters.  Breakup  time 
calculations  are  not  presented  for  droplets  having  initial  diameters  of  10  g;  however,  it  can 
be  inferred  by  linear  extrapolation  of  the  peak  Weber  numbers  (see  Fig.  7)  that  droplets  with 
initial  diameters  greater  than  10  n  and  somewhat  less  than  15  n  would  probably  have 
survived  the  breakup  process  and,  hence,  would  produce  droplets  smaller  than  8.0  in  the 
DET  test  section. 

Further  calculations  could  be  made  to  determine  the  precise  initial  diameter  (between  10 
and  15  /i)  of  the  largest  droplet  that  could  survive  the  breakup  process  described  in  Section 
2.2,  if  such  droplets  were  of  anything  but  academic  interest .  A  more  profitable  pursuit  is 
suggested  by  the  Weber  number  plot  of  Fig.  7,  which  shows  that  the  droplet  breakup 
conditions  occur  only  in  the  nozzle  throat.  Thus,  droplet  breakup  might  be  prevented  if  the 
throat  is  recontoured  to  limit  droplet  acceleration  to  some  low  value  (i.e.,  to  a  Weber 
number  less  than  8.0). 


Nozzle  Station,  in.  Nozzle  Station  in 

Figure  7.  Predicted  droplet  behavior  in  DET  nozzle  throat. 


24 


Normalized  Droplet  Travel  Time, 


A  EDC-TR-79-97 


Nozzle  Station,  in 


Figure  8.  Effect  of  droplet  mass  flow  rate  and  drag 
coefficient  increase  due  to  droplet 
distortion  on  predicted  breakup  time 
for  the  DET  nozzle  throat  region. 
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4.0  DROPLET  ACCELERATION  NOZZLE  CONTOUR  OPTIMIZATION 

The  object  of  the  optimization  process  discussed  here  is  to  limit  the  droplet  Weber 
number  to  a  value  below  which  droplet  breakup  will  not  occur.  A  value  of  8  has  been  chosen 
tentatively  for  this  study;  other  purposes  would  have  produced  other  “optimum”  results. 
Narkis  and  Gal-Or  (Ref.  27)  defined  the  contour  necessary  to  limit  the  evaporation  of  drops 
in  two-phase  flow  and,  with  a  similar  analysis,  came  to  some  strikingly  different 
conclusions,  which  were  correct  for  their  purposes.  In  their  complete  analysis,  Elliot  and 
Weinberg  (Ref.  15)  defined  a  nozzle  contour  necessary  to  maximize  the  exit  velocity  of 
constant-diameter  water  drops.  Other  analyses  of  optimization  criteria  for  two-phase  flow 
systems  may  be  found  in  Refs.  15,  27,  and  28. 

The  salent  feature  of  the  present  analysis  is  that  it  entails  a  simple  modification  of  an 
existing  accurate  predictive  code  that  treats  the  fully  coupled  conservative  equations  for  two- 
phase  flow  of  an  equilibrium  real-gas  and  water  vapor  at  high  pressures  and  moderate 
enthalpy,  including  wall  friction  and  heat  transfer. 

4.1  OPTIMIZATION  CRITERION 


The  basic  optimization  criterion  used  in  this  study— limiting  the  Weber  number  to  an 
arbitrary  value — follows  directly  from  the  conclusions  reached  in  Section  3.0:  that  low 
values  of  Weber  number  associated  with  small,  initial,  drop  diameters  reduce  the  final 
droplet  size  through  increased  evaporation  and  that  high  values  of  Weber  number  result  in 
droplet  breakup.  The  optimization  process  is  involved  at  every  step  in  the  marching  solution 
of  the  system  of  ordinary  differential  equations  of  Section  3.1;  the  process  is  used  to 
determine  the  downstream  pressure  required  to  produce  a  Weber  number  of  the  desired 
limiting  value. 


With  the  assumptions  that  constant  total  enthalpy  and  total  pressure  exist  in  the  gas 
phase  and  that  constant  surface  tension  exists  for  the  predicted  step,  the  downstream 
pressure  is  calculated  from  the  following  functional  dependence: 


/  <30.  <?v, 

p2  =  f(p.Dx,We,Wel,M[.uJ.h1.p1,V  Dr  Dd 

V  1  1  C7X  t?x 


(50) 


This  functional  dependence  [Eq.  (50)]  is  an  approximation,  but  it  works  extremely  well  in 
regions  of  low  pressure  gradient.  The  resulting  value  of  Weber  number,  in  areas  of  high 
pressure  gradient,  is  as  much  as  2,5  percent  low,  which  is  consistent  with  the  approximate 
value  of  the  critical  Weber  number.  Two  more  fundamental  problems  occur  at  the  nozzle 
inlet  and  exit.  At  the  nozzle  inlet  or,  more  precisely,  at  the  droplet  injection  plane,  the 
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optimization  criterion  demands  an  immediate  and  unrealistic  change  in  nozzle  area.  This  is 
circumvented  by  using  the  droplet  injection  conditions  and  the  critical  Weber  number  to 
determine  the  starting  conditions  for  the  optimization  process.  A  similar  problem  occurs 
near  the  exit  of  high  Mach  number  nozzles,  where  the  gas  density  falls  so  low  that  the  local 
Weber  number  is  much  less  than  the  critical  Weber  number.  In  these  cases,  the  optimization 
process  is  no  longer  necessary,  and  the  nozzle  contour  should  be  specified  by  some  other 
criterion.  For  the  conditions  of  this  study,  it  was  found  that  the  droplet  velocity  cannot  be 
appreciably  increased  past  the  point  at  which  the  droplet  critical  Weber  number  is  not  a 
problem. 

4.2  DROPLET  ACCELERATION  NOZZLE 

PERFORMANCE  PREDICTIONS 

A  limited  parametric  study  of  the  factors  affecting  droplet  acceleration  nozzle 
performance  was  made  for  a  HEAT  facility  and  a  few  limited  calculations  were  made  for  the 
DET.  The  HEAT  facility  nozzle  received  more  attention  since  the  characteristically  lower 
expansion  ratio  in  that  facility  allows  faster,  more  economical  calculations  that  retain  all  of 
the  important  characteristics  of  the  processes  common  to  both  facilities,  except  for  droplet 
freezing  that  can  be  predicted  only  by  DET  solutions. 

Preliminary  calculations  indicated  that  droplets  injected  at  velocities  compatible  with  the 
existing  particle  injection  systems  in  the  DET  and  the  HEAT  would  probably  not  survive  to 
the  exit  of  the  optimized  nozzles  at  even  moderate  reservoir  enthalpy  levels.  This  is  further 
aggravated  by  the  fact  that  droplet  and  gas  velocities  vary  with  the  square  of  the  reservoir 
enthalpy. 

One  possible  solution  to  this  problem  is  to  increase  the  droplet  injection  velocity  to 
minimize  the  droplet  mass  loss  in  the  subsonic  section  of  the  nozzle.  As  a  practical  concern, 
the  maximum  droplet  injection  velocity  is  severely  limited  by  the  available  droplet  injection 
pressure  as  is  shown  by  Fig.  2.  The  predicted  effect  of  droplet  injection  velocity  is  shown  in 
Fig.  9.  As  expected,  the  exit  droplet  diameter  increases  almost  linearly  with  increase  in 
injection  velocity  from  13  g  at  600  fps  to  30  g  at  3,000  fps  for  an  initial,  droplet  diameter  of 
50  n .  It  is  unfortunate  that  these  larger  drops  require  very'  long  acceleration  times,  which 
involve  long  nozzle  lengths  and  large  reductions  in  exit  bulk  enthalpy  because  of  wall  heat 
transfer.  The  exit  droplet  velocity  is  4,900  fps  for  an  injection  velocity  of  3,000  ft /sec,  which 
represents  an  increase  of  only  1,900  fps.  Increasing  the  injected  drop  size  from  50  to  150  g 
decreases  the  exit  velocity  from  5,700  to  3,200  fps  (as  shown  in  Fig.  10)  for  injection  velocity 
of  600  fps.  Figure  10  also  shows  that  the  predominant  effect  of  high  injection  velocity  is  to 
increase  the  exit  drop  diameter.  The  exit  bulk  enthalpy  again  decreases  as  the  drop  size 
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Droplet  Injection  Velocity,  fps 

Figure  9.  Effect  of  droplet  injection  velocity  on 
HEAT  optimized  droplet  acceleration 
nozzle  performance. 

increases.  The  performance  of  the  droplet  acceleration  nozzle  for  the  150-/x,  3,000-fps  point 
is  poor  since  viturally  half  the  bulk  enthalpy  is  lost  through  the  nozzle  walls.  The  droplet 
breakup  criterion  severely  limits  acceleration  of  large  droplets. 

Performance  could  be  improved  considerably  if  the  droplet  breakup  could  be  inhibited 
in  some  manner,  such  as  by  increasing  the  surface  tension.  This  is  shown  in  Fig.  11,  in 
which  the  droplet  breakup  Weber  number  is  assumed  to  be  as  high  as  24.  The  increased 
breakup  Weber  number  increases  the  droplet  exit  diameter  and  exit  enthalpy.  The 
performance  of  the  DET  nozzle  with  a  5,000-psi  droplet  injection  system  is  shown  in  Fig.  12. 
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Figure  10.  Effect  of  assumed  droplet  injection  velocity 
on  the  design  point  performance  of  droplet 
acceleration  nozzles  optimized  for  the 
HEAT  facility. 

The  5,000-psi  injection  system  produces  theoretical  droplet  injection  velocity  of  770  fps.  The 
predicted  droplet  exit  velocity  is  rather  insensitive  to  droplet  injection  diameter,  varying  only 
slightly  between  4,700  and  4,800  fps. 


Further  study  is  needed  to  determine  the  full  impact  of  varying  the  other  system  variables 
such  as  throat  diameter  and  nozzle  pressure;  however,  one  problem  that  so  far  seems 
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inherent  in  the  acceleration  nozzle  studied  is  the  large  reduction  in  bulk  enthalpy  attributable 
to  wall  heat  transfer  in  the  excessive  nozzle  lengths.  This  results  in  low  particle  exit  velocities 
even  though  the  particles  are  nearly  in  velocity  equilibrium  with  the  gas.  The  droplet 
velocities  given  here  are  based  on  average  gas  properties;  the  centerline  values  will  be 
somewhat  higher. 


Figure  11.  Effect  of  assumed  critical  Weber  number 

on  the  design  point  performance  of  droplet 
acceleration  nozzles  optimized  for  the 
HEAT  facility. 
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Figure  12.  Droplet  acceleration  performance  predictions 
for  modified  DET  nozzle. 


4.3  DROPLET  INJECTION 

The  tentative  conclusion  that  can  be  made  is  that  high  droplet  injection  velocities  are 
desirable.  A  droplet  injection  velocity  of  3,000  fps  requires  an  injection  pressure  of  60,000 
psi.  Various  investigators  have  used  pressure  intensifiers  to  produce  high  liquid  pressures, 
using  a  low-pressure  gas  source.  Figure  13  is  a  sketch  of  such  a  device  and  Fig.  14  illustrates 
its  application  to  the  HEAT  facility.  The  ratio  of  diameters  of  the  intensifier  can  be  varied  to 
produce  pressure  ratios  of  10:1  or  greater.  Installing  the  intensifier  and  the  arc  heater 
adapter  required  to  minimize  the  friction  losses  will  not  be  a  routine  task. 
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Figure  13.  Droplet  injection  pressure  intensifier. 


Arc-Heated 
Air  Inlet 


Figure  14.  Possible  droplet  injector/arc  heater/acceleration 
nozzle  configuration. 
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5.0  CONCLUDING  REMARKS 

1.  Although  the  analyses  of  the  earlier  test  data  were  not  conclusive,  droplet/ice 
crystal  survival  in  the  DET  can  most  probably  be  attributed  to  high  deformed- 
droplet  drag  loads  and  the  resulting  high  acceleration,  which  reduced  the  droplet 
slip  velocity  before  completion  of  the  breakup  process. 

2.  The  method  of  analysis  described  in  Section  3.1  models  all  of  the  important  flow 
phenomena  in  a  high-temperature  droplet  acceleration  nozzle.  The  water/ice 
phase  transition  should  be  included  in  the  analysis  of  low-temperature  droplet 
acceleration  nozzles  such  as  the  DET. 

3.  Droplet  acceleration  nozzle  contours  can  be  derived  to  inhibit  droplet  breakup 
caused  by  high  Weber  number,  but  these  contoured  nozzles  will  be  very  long  and 
will  have  high  heat -transfer  losses  to  the  wall. 

4.  An  erosion  field  of  35-^-diameter  droplets  can  probably  be  achieved  in  the  DET 
if  a  new  droplet  injection  chamber  and  nozzle  throat  are  designed. 

5.  A  droplet  erosion  field  of  150-^-diameter  droplets  can  probably  be  produced  in 
the  HEAT  facility  if  a  new  droplet  acceleration  nozzle  throat  and  high-velocity 
droplet  injection  mechanism  are  developed. 

6.  Further  analytical  and  experimental  w'ork  should  be  done  concerning  high- 
velocity  droplet  injection  techniques. 

7.  The  one-dimensional  droplet  acceleration  nozzle  code  should  be  modified  to 
include  a  better  wall-heat-transfer  model  and  the  water-ice  phase  change  model. 

8.  A  two-dimensional  droplet  acceleration  nozzle  code  should  be  developed  similar 
to  the  solid,  particle  acceleration  nozzle  code  of  Ref.  28  to  predict  the 
distribution  of  droplets  across  the  nozzle  exit. 

9.  A  full  parametric  study  of  the  droplet  acceleration  nozzle  parameters  should  be 
made  before  any  future  development. 
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APPENDIX 

DERIVATION  OF  TWO-PHASE  FLOW  CONSERV  ATION  EQUATIONS 


The  equations  of  conservation  of  mass,  momentum,  and  energy  used  in  this  study  are 
derived  below  for  one-dimensional  two-phase  compressible  flow  of  a  real  equilibrium  gas 
and  evaporating  liquid  droplets  including  the  effects  of  momentum  and  energy  exchange 
between  the  gaseous  phase  and  the  wall  and  the  gaseous  phase  and  the  droplets. 


The  global  continuity  equation  is  obtained  by  noting 
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From  the  basic  continuity  relationship  mg  =  quA,  it  follows  that 

m  (A- 1 ) 

R  \pdx  u  dx  A  dx  /  dx 

w-here  A  is  taken  to  be  the  total  nozzle  flow  area,  thus  assuming  that  the  droplets  have  no 
volume.  By  defining  a  real-gas  state  parameter  for  the  mixture,  R,  as 


and  noting  that 
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expression  (A-l)  becomes 


l  du  1  dh  1  JR  1  dp  2  dDfj  md  3  dDj 

u  dx  h  dx  R  dx  p  dx  Dx 


m  D,  d  > 

«  d 


(A-2) 


The  momentum  equation  follows  from  Newton’s  Second  Law;  i.e.,  the  change  in 
momentum  is  equal  to  the  sum  of  the  external  forces  applied  to  an  infinitesimal  element: 


since 


dF  =  d(mgu  +  md  Vd) 
dF  =  mgdu  +  (^d  ~u)  d™d  +  ™ 


dV 


dm  =  —  dm  . 
e  d 


since 


then 


dF  =  A  dp  +  jrD^r  dx 

r"gdu  +  (Vd  ~  u)d™d  ~  ™d  dVd  =  AdP  '  n’r)N  rw  d 

du  (Vd_u)  d™J  mJ  dV 


dx 


'd  M1d  u  '  d 

Jx  m  dx 


A  dp  rw 

tn  dx  m 
e  c 


since 


*  Df 


mg 


pu 


r  =  C, 


(A-3) 


then 


1  dmd  _  3  dDd 

m,  dx  D  ,  dx 
o  <1 


du  A  dp  "d  dVd  _  3_^d  <u  -  V  dDd  _  _  2u  Cf 
m  D,  dx  D., 


dx  m  dx  m  d  x 

2  2 


(A-4) 
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The  energy  equation  results  from  the  simple  observation  that  for  a  flowing  system 


d 

dx 


dQ 

W 

dx 


Performing  the  differentiation  and  noting  that 


r*  -DCh  #>u<hw-h0) 

dx 

it  proceeds  directly  that 


d  u  f 

u - l 

fc  T ,  + 

dx 

S  Pd  d 

2 

2  f  m  m  , 

B  d 

nid 

^  +  T  • 

dC 

pd 

+v  dVd  - : 

m  Pd 

~  +  ld 

dx 

dx 

+  d  dx  “ 

Assuming 


A  c 


dx 


and  applying  expression  (A-3)  it  follows  that 


du 

dh  md  / 

'  dVd 

u  —  -f 

—  +  —  l 

^  A  - 

dx 

dx  m  ' 

dx 

dT ,  x  m  ,  ,  dDj 

C  — -  )  +  3  ^  — - - 

f'ddx'  m  D,  dx 


2 


) 


dn 


(h  h  ) 


(A-5) 
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NOMENCLATURE 

A  Cross-sectional  area 

a  Gas  stream  acceleration 

Ac  Gas  stream  acceleration  to  surface  tension  ratio 

Bd  Bond  number 

Cd  Droplet  drag  coefficient 

Cf  Shear  stress  coefficient 

Ch  Stanton  number 

Cp  Specific  heat 

Cq  Flow  discharge  coefficient 

D  Diameter 

Ab  Molecular  diffusion  coefficient 

Dx  Increment  in  x 

F  Force 

g  Gravitational  constant 

h  Specific  enthalpy 

hf  Film  heat-transfer  coefficient 

k  Thermal  conductivity 

Kg  Mass-transfer  coefficient 

M  Mach  number 

m  Mass 

m  Mass  flow'  rate 

Mv  Molecular  weight  of  vapor 

Nu  Nusselt  number  for  heat  transfer 
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Nu' 

P 

Pr 

Q 

q 

q 

R 

r 

Re 

Ru 

5 
Sc 
T 
t 

u 

V 

Vol 

We 

X 

7 

6 
e 
\ 
g 


Nusselt  number  for  mass  transfer 

Pressure 

Prandtl  number 

Normalized  droplet  surface  dynamic  pressure,  total  heat-transfer  rate 

Dynamic  pressure 

Heat-transfer  rate  per  unit  area 

Real-gas  parameter  =  P/gh 

Radius 

Reynolds  number 
Universal  gas  constant 
Surface  area 

Schmidt  number  =  ^/Dabg 
Normalized  droplet  time,  temperature 
Time 

Gas  velocity 
Droplet  velocity 
Volume 
Weber  number 
Nozzle  length  coordinate 
Ratio  of  specific  heats 
Delta  or  differential 
Particle  volume  fraction 


Q 


Latent  heat  of  vaporization 
Micron,  micrometer,  viscosity 
Specific  mass  density 
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a 

Surface  tension 

7 

Shear 

0 

Wall  angle 

SUBSCRIPTS 

am 

Apparent  mass 

B 

Bulk 

b 

Breakup 

Ba 

Bassett 

body 

Denotes  body 

d 

Drop 

g 

Gas 

gr 

Gradient 

i 

Initial 

1 

Liquid 

m 

Average  or  mean 

Max 

Maximum 

N 

Nozzle 

nc 

Noncondensable 

0 

Total 

P 

Surface  pressure 

OR 

Orifice 

r 

Relative 

S 

Stokes 

s 

Surface 
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t  Tip 

v  Vapor 

w  Wall 

oo  Infinity 
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